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We present lattice results for the form factors relevant in the K — * 7t£vi and D — > nlvt semileptonic 
decays, obtained from simulations with two flavors of dynamical twisted-mass fermions and pion 
masses as light as 260 MeV. For K — > nlv decays we discuss the estimates of the main sources 
of systematic uncertainties, including the quenching of the strange quark, leading to our final re- 
sult /+(0) = 0.9560(57) s tat.(62)syst.- Combined with the latest experimental data, our value of 
/+(0) implies for the CKM matrix element |V MJ | the value 0.2267 (5)exp. (20) ^ /q\ consistent with 
the first-row CKM unitarity. For D — > %lVt decays the application of Heavy Meson Chiral Per- 
turbation Theory allows to extrapolate our results for both the scalar and the vector form factors 
at the physical point with quite good accuracy, obtaining a nice agreement with the experimental 
data. In particular at zero-momentum transfer we obtain /+(0) =0.64(5). A preliminary analysis 
of the discretization effects is presented and discussed. 



The XXVII International Symposium on Lattice Field Theory - LAT2009 
July 26-31 2009 

Peking University, Beijing, China 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Semileptonic Form Factors from Two-Flavor Lattice QCD 



S. Simula 



1. Introduction 



Semileptonic decays of pseudoscalar mesons can provide important information on the weak 
mixing of quark flavors, which in the Standard Model constructs the well-known CKM matrix 
[j|]. In order to extract from the experimental data precise values of the relevant CKM entries it is 
necessary to determine precisely the matrix elements of the weak hadronic current. 

In the case of the semileptonic decay H — ► Pl\g the matrix element of the weak vector current 
can be written in terms of two form factors, the vector, f+(q 2 ), and the scalar, fo{q 2 ) ones, namely 

(P(p P )\V^\H(p H )) = (p P + p H -Aff + (q 2 )+A^fo(q 2 ), (1.1) 

where A = q (Mf\ — Mp) / q 2 and q = ph — Pp is the 4-momentum transfer. 

In this contribution we present the lattice results for the vector and scalar form factors obtained 
by the European Twisted Mass (ETM) Collaboration in the case of the K — ► 7l£vi and D — > n£ve 
semileptonic decays, which are relevant for the determination of the CKM matrix elements \V US \ 
(known also as the Cabibbo's angle) and \V cc i\, respectively . 



2. K — > %£v decays 

The relevant hadronic quantity in the case of the K — ► %£v decays is the vector form factor at 
zero-momentum transfer, /+(0). Its first determination dates back to the eighties, i.e. to the work 
of Ref. in which Chiral Perturbation Theory (ChPT) and the quark model were employed. 

The determination of / + (0) using lattice QCD started only more recently with the quenched 
calculation of Ref. [Qj, where it was shown how /+(0) can be determined at the physical point 
with ~ 1% accuracy. The findings of Ref. [||] triggered various unquenched calculations of /+(0), 
namely those of Refs. [Q, ||, ^] with Nf = 2 and pion masses above ~ 500 MeV and the recent one 
of Ref. [|7|] with Nf = 2 + 1 and pion masses starting from 330 MeV. 

In Ref. [JSj] a new lattice result for /+(0), namely 

/+(0) = 0.9560 ±0.0057 stat . ±0.0062 syst . = 0.9560 ±0.0084 , (2.1) 

was obtained by the ETM Collaboration using gauge configurations with Nf = 2 flavors of dynam- 
ical twisted-mass quarks [^] and simulating pion masses from 260 MeV up to 575 MeV. 



Our new determination ( |2.1| ) agrees very well with the Leutwyler-Roos result [|2p and with 
previous lattice calculations at Nf = [§], Nf = 2 [@, |, g| and Nf = 2 + 1 [@|. Using the latest 
experimental determination of the product |V r MiS |/ + (0) = 0.21668(45) [|l^, [IT]] we get from (2.1) 



|K.s-| = 0.2267 ± 0.0005 eX p. ± 0.0020 /+(0) . (2.2) 



Combining this value with \V ud \ = 0.97418(27) and \V ub \ = 0.00393(36) from PDG2008 ^ the 
first-row CKM unitarity relation becomes 

\V ud \ 2 + \V us \ 2 + \V ub \ 2 = 1.0004 ±0.0015. (2.3) 

Our final value ( pT| ) includes the estimates of all sources of systematic errors: discretization, 
finite size effects (FSE's), q -dependence, chiral extrapolation and the effects of quenching the 
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strange quark. In Ref . [fSf] the chiral extrapolation and the related uncertainty on / + (0) were inves- 



tigated using both SU(3) and, for the first time, SU(2) ChPT Q12|], obtaining fully consistent results. 
Also the q 2 -dependence of the form factors was investigated by considering different functional 
forms. The systematic error associated both to the g 2 -dependence and to the chiral extrapolation 
was determined quite accurately and turned out to be 0.0035 [Jsj] . 

We illustrate now in more details the estimates of the remaining sources of systematic errors, 
namely finite size, discretization and the quenching of the strange quark. 

Finite Size. We have performed simulations close to M n ~ 300 MeV using two lattice volumes, 
24 3 • 48 a 4 and 32 3 • 64 a 4 , for a lattice spacing equal to a ~ 0.088 fm. The two simulations 
correspond to M n L ~ 3.2 and 4.2, respectively. As described in Ref. ||8|], a smooth interpolation of 
/+ (0) at the physical strange quark mass can be obtained by fixing the combination (2M§ — M 2 ) at 
its physical value. Thus for each pion mass M n a reference kaon mass Mg is defined as 



2[m; 



K f } 2 



■Mi 



2[M p K hys ] 2 



[Ml 



physYl 



(2.4) 



with Ml 



phys 



135.0 MeV and M 



phys 
K 



494.4 MeV. 



The results for /+(0), obtained adopting either the pole-dominance or a quadratic fit for de- 
scribing the ^-dependence of the form factors (see Ref. [^]), are shown in Fig. [j](a) versus the 
lattice size L/a. For matrix elements like (ti\V^\K), involving one particle in the final states, FSE 
are known to be exponentially suppressed. Assuming a dependence of the form A + Be~ M " L /L 3 / 2 
the residual FSE, corresponding to the difference between the value at infinite volume and the one 
calculated at the largest lattice volume, turns out to be equal to 0.0018. 



1.000 



0.995 v 



0.990 



0.985 



0.980 



0.975 



a ~ 0.088 fm 




1.0000 



T— I— I— I— |— 1— I— I— I— J— I— I— I— I— |— I— I— I— 

(a) M ~ 300 MeV, M = M ref ~ 530 MeV '. 

v ' it K K 

• pole fit 
□ quadratic fit 

A + B exp(-M L) / L 3/2 j° 0.9950 



0.9975 



0.99Z5 



0.9900 




• pole fit 
□ quadratic fit 

linear fit in a 2 



_l_ 



_l_ 



20 



25 



30 

L/a 



35 



40 



0.000 



0.004 



0.008 



0.012 



a 2 (fm 2 ) 



Figure 1: Vector form factor at zero-momentum transfer, /+(0), versus the lattice size in lattice units (a) and 
the squared lattice spacing (b). The values of the pion mass are reported in the inset, while the kaon mass 



is fixed at the corresponding reference values given by Eq. (2.4). In (a) the dotted line represents the value 
/+(0) — A in the limit of infinite volume. 



Discretization. We have performed simulations dXM n ~ 470 MeV using three lattice spacings: 
a ~ 0.069,0.088 and 0.103 fm. The results for /+(0), shown in Fig. [Jb), exhibits a clear, linear 
(in a 2 ) increase toward the continuum limit, consistent with the automatic 0{a) -improvement at 
maximal twist [jl3p. The difference between the value in the continuum limit and the one at a ~ 
0.088 fm is equal to 0.0037, which represents our estimate of the contribution of discretization 
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effects to the systematic error in Eq. (2.1). A complete study of the scaling property of /+(0) at 
various pion masses is in progress. It will allow us to compute the continuum limit, reducing in 
this way significantly the error due to lattice artifacts. 

Quenching of the strange quark. The effect of our partially quenched (PQ) setup can be esti- 
mated within SU(3) ChPT, which provides a systematic expansion of /+(0) of the type 



/+(0) = l+/2 + /4 + /6 + . 



(2.5) 



where = &[Mj, J (4nf n ) n ] and the first term is equal to unity due to the vector current conser- 
vation in the SU(3) limit. 

Because of the Ademollo-Gatto theorem Q], valid also in both quenched (Q) [||] and PQ JT5| ] 
setups, the first correction fa does not receive contributions from the local operators of the effective 
theory and can be computed unambiguously (for any Nf) in terms of the kaon and pion masses and 
the pion decay constant f n . At the physical point it takes the values: f% = +0.022 in the quenched 
-0.0168 for our PQ setup with N f 



case N 



2 [15] and f 2 



0.0226 for N f = 2 + 1 



]. Thus the effect of quenching the strange quark is exactly known at NLO: fa — p 



PQ 



-0.0058 



(~ 26% of fx). This correction was taken into account in Ref. [g] and it has no error. Note that the 
difference between the values of fa at Nf ■ = 2 + 1 and Nf = 2 is almost an order of magnitude less 
than the difference between those at Nt = 2 + 1 and Nf = 0. In our opinion this should be traced 
back to the facts that fa is dominated by meson loops and the pion contribution is the same in the 
Nf = 2 and Nf = 2 + 1 theories. 

The task is thus reduced to the problem of estimating the quenching effect on the quantity 



A/ = / 4 + / 6 + ... = / + (0)-(l+/ 2 ) 



(2.6) 



The results obtained for Ap by the ETM Collaboration at N f = 2 [@| and by the RBC/UKQCD one 
at Nf ■ = 2 + 1 |^] are compared in Fig. ||. It can clearly be seen that the effect of quenching the 
strange quark is well within the statistical uncertainties found by the two Collaborations. 
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Figure 2: Values of the &(p 6 ) term Af [Eq. (fj)/ obtained by the ETM /§/ and RBC/UKQCD /@/ Collab- 
orations taking into account the values of the NLO term fa appropriate for N f = 2 and Nf = 2 + 1. 



In Ref. [Q] the relative quenching error on Af has been estimated to be at most 50% of the 
same relative effect on fi. Such an estimate is based on the observation that, while the NLO term 
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fi is expected to be sensitive to the number of sea-quark flavors being only determined by the 
contribution of meson loops, the &{p ) term Af receives important contributions from the local 
terms of the effective theory, which are expected to be dominated by the physics of the nearest 
resonances. Our estimate corresponds to a systematic error of 0.0028 (i.e., ~ 13% of Af), which 
incidentally turns out to be of the same size of the difference between the ETM result for Af at 
Nf = 2 and the quenched one of Ref. [[|]. Thus we expect our estimate of the quenching error to be 
a quite conservative one. 

To close this Section we have collected the budget for the systematic error of the ETM result 
in Table |. 



Source 


systematic error 


%of[l-f + (0)] 


q 2 — dependence and chiral extrapolation 


0.0035 


8 


finite size 


0.0018 


4 


discretization 


0.0037 


8 


quenching of the strange quark 


0.0028 


6 


Total {in quadrature) 


0.0062 


14 



Table 1: Budget of the systematic error for the ETM determination f |2.7| ) of the vector form factor at zero- 
momentum transfer, /+ (0), obtained in Ref. /]s]/. 



3. D — ► n£v( decays 

In the case of the semileptonic decays of a heavy meson H it is convenient to use a decompo- 
sition of the matrix element of the weak vector current in which the form factors are independent 
of the heavy-meson mass M# in the static limit, namely 



(P(p P )\V^\H(p H )) = ^2M H [v" f v (E) + p^f p (E)] , (3.1) 

where v = ph/Mh, Px = pp—Ev and E = v • pp = {Mjj +Mp — q 2 ) /2Mh is the energy of the final 
meson in the rest-frame of the initial one. The relation between the form factors f v (E) and f p (E) 



and those appearing in Eq. (1.1) is 



f+(q 2 ) = [ME) + (Mh — E)f p (E)] /^/m; , (3.2) 
fo(q 2 ) = [(M H -E)f v {E) + {E 2 -M 2 )f p {E)} ^2M H '/ \M 2 - M 2 ) . (3.3) 

In the static limit both the mass and the energy dependence of the form factors f v (E) and 
f p (E) have been investigated within the Heavy Meson ChPT (HMChPT) in Ref. [16]. For a pion 
in the final state, i.e. P = n, one has at NLO 

ME) = Do [1 +D l {E)M 2 n + D 2 (E)- 3(1 + 3g 2 )M 2 K L(M 2 K )/4 

- 2{E 2 -M 2 n ) L(Ml) - 2M n E F(E/M n )] , (3.4) 
f p (E) = C [l+C 1 (E)Ml + C 2 (E)-3(\ + 3g 2 )MlL(Ml)/4\ /(E + A*) , (3.5) 
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where Q and D, (i = 0, 1,2) are unknown low-energy constants (LEC's), g is the H*Hn coupling 
constant (with H* being the vector resonance of the heavy meson H), A* = Mr* —Mh, L(M 2 ) = 
log(Ml)/(4Kf n ) 2 and F(x) = 2^/x 2 - 1 log[x+ y/x 2 - \}/{Anf n ) 2 for x > 1. 

In addition a generalization of the Callan-Treiman relation [|l7]] constrains the value of fo{q 2 nax ) 
to be equal to the ratio of the leptonic decay constants fjj / f n at the chiral point M n = 0. Therefore, 
using SU(2) HMChPT for f H and SU(2) ChPT for f K , one gets at NLO 



fH/fx = \/VM^D [l+D 2 (p)] [\+BM 2 K + {5-9 g 2 )M 2 K L{Ml)/A\ , (3.6) 

where B is an unknown LEC. 

We have calculated the vector and scalar form factors f+(q 2 ) and fo(q 2 ) as well as the decay 
constants fy and f n JT8|], using the gauge configurations generated by the ETM Collaboration [Q] 
with Nf = 2 flavors of dynamical twisted-mass quarks at a single lattice spacing a ~ 0.088 fm for 
various values of the sea quark mass. The valence light-quark mass is kept equal to the sea quark 
mass to get unitary pions with a simulated mass ranging from ~ 260 to ~ 575 MeV, as in the study 
of K — > n£vc decays. For each pion mass we use three values of the charm quark mass to allow for 
a smooth, local interpolation of our results to the physical D-meson mass. At the two lowest pion 
masses the lattice volume is L? ■ T = 32 3 • 64 a 4 , while at the higher ones it is 24 3 • 48 a 4 in order to 
guarantee that M n L > 3.7. 



We have then applied Eqs. (3.5) and (3.6) for a simultaneous fit of the energy and pion mass 



dependence of our results. Simple polynomial parameterizations of the energy dependence of the 



LEC's C12 and D\ 2 have been adopted, while the value g = 0.6 has been taken from Ref. [19] 
and the quantity A* has been fixed at its value at the physical point (A* = 138 MeV) . The range 
of values of q 2 covered by our data is quite large, extending from q 2 up to q 2 = q 2 max , which 



corresponds to values of the energy E up to « 1 GeV. The quality of the fit provided by Eqs. ( p.5| ) 



turns out to be remarkably good, though the chiral expansion of Ref. [ |16| ] is in principle limited to 
the static limit and to values of the energy E well below the scale of chiral-symmetry breaking. 

The extrapolation of the vector and scalar form factors to the physical pion mass is quite 
accurate in the full range 

tfmax> as shown in Fig. ^. Our results are also in good agreement 
with the latest experimental data from the CLEO Collaboration |p0[|, obtained after assuming for 
the CKM matrix element \V c d\ the value implied by unitarity. Only around q 2 « q 2 max our lattice 
predictions for/ + (g 2 ) are slightly below the experimental results. Atg 2 = we get /+ (0) =0.64(5) 



which agrees with the lattice result / + (0) = 0.64(3) (6) obtained in Ref. Q with N f = 2 + 1. 

Since our results have been obtained at a single value of the lattice spacing (a ~ 0.088 fm) and 
a heavy mass, like the charm one, is involved, the question of possible sizable discretization effects 
naturally arises. We have therefore computed the form factors for M n ~ 470 MeV at three values 
of the lattice spacing (a ~ 0.069,0.088 and 0.103 fm). The results are shown in Fig. |j. 

It can be seen that discretization effects are small around ^0 (at the level of the percent 
between the two finest lattices), while they increase toward q 2 = q 2 ^, particularly in the case of 
the vector form factor f+(q 2 ). Thus we observe that discretization effects related to the presence of 
the charm quark mass are quite limited for our setup. This may be related to the fact that the form 
factors are extracted from ratios of correlation functions, in which lattice artifacts may partially 
cancel out. We observe moreover that the discretization effects on /o(^mar) we °f tne order of few 
percent, i.e. similar to those found for fo/fn in the case of our action JT~8|] . 
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Figure 3: Vector [f+(q 2 ]] and scalar [fo(q 2 )] form factors for the D — > 7t£ve decay versus the squared 
4-momentum transfer q 2 . The bands correspond to the regions selected at 1(7 level by the chiral fit (3.5) 
applied to our lattice results. The dots and the squares are the experimental data from Ref. /p0|/. 
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Figure 4: Vector (blue markers) and scalar (red markers) form factors for the D — > 7t£"Vi decay versus the 
squared 4-momentum transfer q 2 in units of the Sommer parameter ro, obtained at three lattice spacings for 
M n ~ 470 MeV and smoothly interpolated at the physical D-meson mass. 



A complete study of the scaling property of the vector and scalar form factors at various 
pion masses is in progress. The results presented in Fig. |], however, makes us confident that 
the agreement with the experimental data visible in Fig. || will not be spoiled by a more detailed 
analysis of discretization effects. 
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